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ABSTRACT 


We have conducted the first parallax and proper motion measurements of 
6.7 GHz methanol maser emission using the Australian Long Baseline Array 
(LBA). The parallax of G 339.884—1.259 measured from five epochs of observa¬ 
tions is 0.48±0.08 mas, corresponding to a distance of 2.11 q 3 kpc, placing it in 
the Scutum spiral arm. This is consistent (within the combined uncertainty) with 
the kinematic distance estimate for this source at 2.5±0.5 kpc using the latest 
Solar and Galactic rotation parameters. We find from the Lyman continuum 
photon flux that the embedded core of the young star is of spectral type Bl, 
demonstrating that luminous 6.7 GHz methanol masers can be associated with 
high-mass stars towards the lower end of the mass range. 

Subject headings: astrometry - Galaxy:kinematics and dynamics - masers - 
stars: formation 



1. Introduction 


To properly understand the Milky Way’s scale and shape as well as the physical 
properties of the objects within it, including their mass, luminosity, ages and orbits, 
we need to be able to accurately measure distances to Galactic sources. Unfortunately, 
this fundamental attribute is one of the most difficult measurements to make. Common 
techniques include the use of a Galactic rotation model, or “standard candles” where 
the distances to objects with known luminosities can be determined from their observed 
brightness. These indirect methods of measurement can result in large uncertainties unless 
an accurate, reliable and robust model is available. For example Xu et al. (2006) found that 
the kinematic distance was factor of 2 greater than the parallax distances to W3(OH) due 
to the peculiar velocity of the source. Uncertainties and errors in distance determination 
naturally propagate into the estimation of other important properties, such as luminosity 
(L oc D 2 ) and mass (L oc M a , where 3 < a < 4 for main sequence stars), and hence there 
is a need to constrain any errors in distance as much as possible. One of the most direct 
methods to determine distances beyond our solar system is through the use of trigonometric 
parallax. 

A decade ago, Honma et al. (2004) used the VLB1 Exploration of Radio Astrom¬ 
etry (VERA) array to conduct high precision astrometry on the water maser pair 
W49N/OH43.8—0.1 and achieved an accuracy of 0.2 mas. This success demonstrated the 
feasibility of performing maser astrometry on Galactic scales and subsequently several 
groups obtained parallax results to masers and their associated star formation regions 
(Hachisuka et al. 2006; Xu et al. 2006; Honma et al. 2007). 

Following this, the BeSSeL project was launched as a comprehensive Northern 
hemisphere survey to measure accurate distances to high-mass star formation regions 
(HMSFR) and associated Hll regions of our Milky Way galaxy using trigonometric parallax 
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(Reid et al. 2009a; Moscadelli et al. 2009; Xu et al. 2009; Zhang et al. 2009; Brunthaler et al. 
2009; Reid et al. 2009b; Sanna et al. 2009; Brunthaler et al. 2011). This is being undertaken 
using the National Radio Astronomy Observatory’s (NRAO) Very Long Baseline Array 
(VLBA) to measure position changes of methanol and water masers in the Galactic disk, 
with respect to distant background quasars. Phase referencing of one of the quasar and 
maser emission data to the other (Alcf 1988; Beasley & Conway 1995; Reid & Honma 
2014), combined with careful calibration of atmospheric effects, allows the change in the 
relative separation between the masers and quasars to be accurately measured. Repeated 
observations timed to maximise the measured amplitude of the parallax signature can 
measure the parallax to an accuracy of up to 10 /rarcsec and simultaneously determine 
the source proper motions to ~1 /xarcsec year -1 (Reid & Honma 2014). Thus far the 
combination of all astrometric VLB1 observations (including the European VLBI Network 
(EVN))(e.g. Reid & Honma 2014; Chibueze et al. 2014; Imai et al. 2013; Honma et al. 
2012; Sakai et al. 2012; Rygl et al. 2010, 2008) has yielded more than 100 parallax and 
proper motion measurements for star forming regions across large portions of the Milky 
Way visible from the Northern hemisphere. Having determined the position at a reference 
epoch, parallax and proper motion of the masers, the complete three-dimensional location 
and velocity vectors of these sources relative to the Sun can be found. Reid et al. (2014) 
used these measurements to fit a rotating disk model of the Milky Way and from this data 
were able to refine the best-fit Galactic parameters finding the circular rotation speed of 
the Sun, © 0 = 240 ± 8 km s^ 1 and distance to the Galactic centre R 0 = 8.34 ± 0.16 kpc. 

To date, parallax distances to masers have only been measured using Northern 
hemisphere VLBI arrays, and therefore, the sources for which accurate distances have been 
measured are heavily concentrated towards objects in the first and second quadrants of the 
Galaxy. The sources with trigonometric parallax measurements compiled in Reid et al. 
(2014) predominantly lie within the Galactic longitude range of 0° < l < 240°. There are 



two sources within the fourth quadrant of the Galaxy, with the most southerly of these 
at a declination of —39°. Sources at southerly declinations achieve only low elevations 
for Northern hemisphere antennas, and for these, the measurement error in the zenith 
tropospheric delay estimate produces significant degradation in the astrometric accuracy, 
and corresponding parallax measurement (Honma et al. 2008). Therefore, the updated 
estimates for fundamental Galactic parameters such as ©o and i? 0 are based on data which 
are restricted to only a little over half the total Galactic longitude range. In order to ensure 
that models of Galactic structure and rotation are reliable, they must be derived from a 
more uniformly sampled distribution including sources in the third and fourth quadrants. 

Distances to southern maser sources and their associated Hu regions have exclusively 
been determined via indirect methods, such as through kinematic distance estimates (e.g. 
Caswell et al. 1975, 2010), with the ambiguity being resolved with Hi absorption against Hu 
continuum (Jones & Dickey 2012; Jones et al. 2013), Hi Self-absorption (HlSA) (Green & 
McClure-Griffiths 2011) and radio recombination lines (Sewilo et al. 2004). However, there 
are issues which can produce significant errors in these techniques including, potentially 
large peculiar motions, broad spectral line profiles, and the as yet unknown relationship 
between cold Hi distribution against the maser associated component in HlSA emission. 

In 2008 a project was initiated to precisely measure the positions of maser sources 
relative to extragalactic quasars using the Australian Long Baseline Array (LBA). The aim 
of this ongoing project is to determine the parallax distances to 30 prominent southern 
HMSFRs. Current maser parallax measurements made with the VLBA and VERA have 
primarily been towards either 12.2 GHz methanol masers or water masers. However, not 
all LBA antennas have receivers capable of observing at 12.2 GHz, and many of the LBA 
antennas have significantly poorer sensitivity at 22 GHz than they do at frequencies less 
than 10 GHz. Hence the 6.7 GHz class II methanol maser transition was considered to 
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provide the best targets for parallax observation with the LBA. They have strong, stable 
and compact emission over the timescales required for parallax measurements, and are well 
sampled from the Southern hemisphere, with close to 1000 sources being documented in the 
Methanol Multibeam (MMB) survey (Caswell et ah 2010; Green et al. 2010; Caswell et ah 
2011; Green et al. 2012, a sensitive, unbiased search for 6.7 GHz methanol masers in the 
southern Galactic plane). Here we present the first trigonometric parallax measurements of 
a southern 6.7 GHz methanol maser. 


2. Overview of G 339.884-1.259 

G 339.884—1.259 is one of the strongest (1520 Jy at —38.7 km s -1 (Caswell et al. 
2011)) 6.7 GHz methanol masers and has been well studied in a range of maser and thermal 
molecular transitions (e.g. Norris et al. 1993; De Buizer et al. 2002; Ellingsen et al. 2004, 
2011) and in continuum emission at a range of wavelengths (e.g. Ellingsen et al. 1996; Walsh 
et al. 1998; Ellingsen et al. 2005). The autocorrelation spectrum of the 6.7 GHz methanol 
maser emission of this source from the 2013 June epoch of observation is shown in Figure 1. 
Its peak flux density and spectral profile has been relatively stable over the last 20 years, and 
interferometric observations of the masers indicate several compact features (Norris et al. 
1993; Ellingsen et al. 1996; Caswell et al. 2011), making it a suitable candidate for phase 
referenced astrometry. Previous estimates of the distance to G 339.884—1.259 have utilised 
the kinematic distance method (e.g. Caswell & Haynes 1983; Green & McClure-Grifhths 
2011) which yields a kinematic heliocentric distance of between 2.5 to 3.0 kpc (depending 
on the assumed Galactic parameters). The coordinates used for G339.884—1.259 in the 
current observations are given in Section 5.3 with updated coordinates presented in Table 1. 

Methanol masers in G 339.884—1.259 were first observed by Norris et al. (1987) at 
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Table 1: Coordinates of observed sources. The separation and position angle columns de¬ 
scribe the offset between the respective quasar and G 339.884—1.259 in the sky. The reported 
positions of G 339.884—1.259 and J 1654—4812 have been revised (see Section 5.3) based on 
the 2013 March epoch. We failed to detect J 1648—4826, J 1644—4559, J 1648—4521 and 
J 1649—4536 in our observations. The upper limits for detection is 5 times the image RMS 
(from a box of size 1.5 x 1.5 arcsec). 


Source 

Correlated 

Separation 

Position 

RA 

Dec 


flux 


angle 




(mJy) 

(°) 

(°) 

(h m s) 

( ° , H) 

Maser: 






G 339.884-1.259 

— 

— 

— 

16 52 04.6776 

-46 08 34.404 

Detected quasars: 






J 1706-4600 

131.1 

2.48 

88.1 

17 06 22.0503 

-46 00 17.824 

J 1654-4812 

3.8 

2.11 

169.9 

16 54 18.2448 

-48 13 03.756 

Non-detected quasars: 






J 1648-4826 

<0.5 

2.36 

193.3 

16 48 47.9200 

-48 26 18.800 

J 1644-4559 

<0.5 

1.27 

276.5 

16 44 49.2856 

-45 59 09.646 

J 1648-4521 

<0.9 

1.02 

319.1 

16 48 14.2110 

-45 21 38.090 

J 1649-4536 

<0.4 

0.72 

317.1 

16 49 14.7810 

-45 36 31.190 
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12.2 GHz and high resolution synthesis images were made of the emission at 6.7, and 
12.2 GHz by Norris et al. (1993). Interestingly, the maser emission was found to have a 
linear spatial distribution and had a corresponding monotonic velocity gradient. In light of 
this, G 339.884—1.259 became a prime candidate for astronomers to model the conditions 
for high-mass star formation (Ellingsen et al. 1996; Norris et al. 1998; Phillips et al. 1998; 
Walsh et al. 1998; De Buizer et al. 2002; Dodson 2008). Ellingsen et al. (1996) proposed 
that the masers are located within a cicrumstellar disk, and Norris et al. (1998) showed 
that the emission fit a model of a Keplcrian disk around an OB type star. However, these 
claims were disputed by De Buizer et al. (2002) who demonstrated that what was initially 
thought to be a single circumstellar disk could be resolved into three mid-infrared sources 
near the location of the methanol masers. Dodson (2008) attempted to test the hypothesis 
of a circumstellar disk in G 339.884—1.259 by making polarization measurements of the 
magnetic fields associated with this source. From their images, they report a primarily 
disordered field accompanying much of the emission, and proposed that this matches the 
expectations for the masers being associated with an outflow-related shock. One small 
region of methanol maser emission does show a magnetic field direction perpendicular to 
the elongation of the maser emission, suggestive of a disk (Dodson 2008). However, the 
mass of the enclosed object (assuming Keplarian rotation) is of magnitude of only 0.03 M 0 
(assuming a distance of 3 kpc to G 339.884—1.259). 

A large number of additional class II methanol maser transitions have been detected 
towards G 339.884—1.259. Observations of multiple maser transitions are required to place 
constrains on the environmental conditions in the region surrounding the young high-mass 
star where the masers arise. This is because any model will need to account for the specific 
conditions required for the observed transition (Gragg et al. 1992; Sobolev et al. 1997). 
Norris et al. (1993) made high-resolution maps of the 6.7 and 12.2 GHz maser spots in 
G 339.884—1.259 and found little evidence for spatial coincidence with the spots at different 
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frequencies. Caswell et al. (2000) made observations of 107.0 and 156.6 GHz methanol 
masers in G 339.884—1.259 and found the emission peak at these frequencies coincided with 
the 6.7 GHz maser site to within 5". These are extremely rare transitions, with only 22 
and 4 detections respectively from a pool of 80 sources. Ellingsen et al. (2004) discovered 
emission from the 19.9 GHz transition in G 339.884—1.259, and in Krishnan et al. (2013) 
it was shown that there was little evidence for spatial or velocity coincidence between 
the masers at 6.7 and 19.9 GHz, with only 2 maser components identified at 19.9 GHz 
as opposed to ~10 components in 6.7 and 12.2 GHz observations with similar angular 
resolution (Norris et al. 1993). 

Ultracompact (UC) Hll regions are bubbles of ionized gas associated with newly 
formed massive stars, and the first detection of such a region in G 339.884—1.259 was by 
Ellingsen et al. (1996). The emission was measured at 8.5 GHz and had a peak brightness of 
6.1 mJy beam -1 . The peak of the methanol emission at —38.7 km s -1 was found to be offset 
from the continuum peak by 0.6", with the methanol masers lying in a line approximately 
across the centre of the continuum emission and in an orientation which is perpendicular to 
the direction of extension of the UCHlI region. De Buizer et al. (2002) interpret the radio 
continuum emission is being due to an ionized outflow along the axis of extension. 


3. Observations 

The Long Baseline Array (LBA) is a heterogeneous VLBI array with either 5 or 6 
antennas available for the observations reported here. In the period spanning 2012 March 
to 2014 March, a total of six epochs of observations were undertaken (LBA experiment 
code v255, epochs q to v inclusive) towards the southern 6.7 GHz methanol maser 
G339.884—1.259 (see Table 2). The LBA antennas available for one or more epochs were 
the Australia Telescope Compact Array (ATCA), Ceduna 30m, Hartebeesthoek 26m, 
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Hobart 26m, Mopra 22m and Parkes 64m antennas. The ATCA is itself a connected 
6-element interferometer, which was operated in tied-array mode, with the outputs from 
either four or five 22 metre antennas phased and combined. 

The observations for each epoch typically lasted for between 18 to 24 hours, of which 
approximately one-third of the time was utilised for observation of G 339.884—1.259 (and 
associated calibration observations). The remaining time was used for observations of other 
sources, the results of which will be reported in future publications. The setup consisted 
of two different frequency configurations to maximize on the sensitivity requirements of 
the different modes. The first configuration was to record continuum observations for 
calibration of the tropospheric component of the delay, and the second for phase referenced 
maser observations. 

The first configuration utilized as broad a frequency range as possible. The 
heterogeneous nature of the array meant that with different receiver front ends, it was 
necessary to compromise the frequency setup so that it could be adopted at all antennas. 
The LBA Data Acquisition System (DAS) can record the observed signals onto two 
independent IF bands. The optimal frequency configuration which is able to accommodate 
these restrictions is to have 4 x 16 MHz bands, the first pair (LBA DAS IF 1 with RR 
polarization) with lower-band edges at 6300 and 6316 MHz and the second pair (LBA 
DAS IF 2 with LL polarization) with lower-band edges at 6642 and 6658 MHz. Brief 
observations of approximately 12 to 18 quasars (~2 minutes per source) with point-like 
structure at high resolution were undertaken for a broad azimuth range (generally at low 
elevation). These quasars were selected from the International Celestial Reference Frame 
(ICRF) Second Realization catalogue (Ma et al. 2009). These “ICRF observations” were 
used to determine the troposphere path length contributions to the delay (r), as well as to 
model the clock drift rate at the observatories (Section 4.1 and Appendix A). The ICRF 
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observations were grouped into 45 minute blocks with an interval of between 3 to 6 hours 
separating consecutive blocks. 

The second frequency configuration was for phase referencing between the maser 
G 339.884—1.259 and a nearby background quasar for parallax determination. The LBA 
DAS system was set to record dual circular polarisation for 2 x 16 MHz bands with 
lower-band edges at 6642 and 6658 MHz. The phase referencing technique was employed 
by alternating scans for 2 minutes on the target maser with scans lasting 2 minutes 
on a nearby quasar. In order to achieve accurate phase referenced astrometry, suitable 
background quasars had to fulfill the criteria that they have little or no extended structure, 
their coordinates be known to an uncertainty of <1 mas and that they are in close angular 
proximity to the associated maser of ~1° (Reid et al. 2009a). The primary databases 
which were used in our search for quasars were the AT20G (Hancock et al. 2011), Astrogeo 
(Petrov et al. 2011) and ATPMN (McConnell et al. 2012) catalogues. A list of potential 
background quasars which were observed in conjunction with G 339.884—1.259 is given 
in Table 1. During the data reduction process, the phase information from the maser 
emission in a single channel (where the emission is strong and compact) is transferred to 
the nearby quasar (see Section 4.2). In doing this, it is assumed that the state of the array 
has remained constant in the time interval between consecutive maser scans, so that the 
quasar phase can be interpolated with small error contributions (Fomalont 2013). 

The data were correlated with the DiFX correlator at Curtin University (Dcllcr et al. 
2011). As the maser emission covers only a small fraction of the 32 MHz of recorded 
bandwidth, it is not necessary to correlate the entire band. Hence, for the maser data 
a 2 MHz zoom-band was correlated with 2048 channels (1 MHz over 1024 channels for 
the 2014 March epoch), giving a spectral resolution of 0.977 kHz corresponding to a 
velocity separation of 0.055 km s -1 . For the observations of the background quasar the 
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entire observed bandwidth was correlated. In the 2012 March and 2013 March epochs, 
256 spectral channels were used per 16 MHz bandwidth, corresponding to a resolution of 
62.5 kHz. In the remaining epochs, 32 spectral channels were used per 16 MHz bandwidth, 
corresponding to a resolution of 500 kHz. The ICRF data were correlated using the same 
spectral resolution as the background quasar observations in each epoch. 


4. Calibration 

We used the Astronomical Image Processing System (AIPS) (Greisen 2003) for data 
processing, employing the same data reduction steps across all epochs and based on the 
procedure described in Reid et al. (2009a). 


4.1. ICRF Data 

Prior to using observations of the ICRF catalogue quasars to determine the tropospheric 
and clock delay parameters for each antenna, it was necessary to remove the estimated 
ionospheric delay (determined from global models based on GPS total electron content 
(TEC) observations (Walker & Chatterjee 1999)), the Earth Orientation Parameters 
(EOPs) and parallactic angle effects. We then performed delay calibration, using a strong 
ICRF source with well known position (<1 mas) (Ma et al. 2009) (see Table 3 for details). 
A least squares fit based on the approach in Reid et al. (2009a) was then employed to 
determine the zenith atmospheric delay and the results (detailed in Appendix A) were 
applied to the phase referenced observations in Section 4.2. 
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Table 2: Stations which participated in the observations of G 339.884—1.259 and 


J 1706-4600. 


D.O.Y 

Start UT 

Code 

Epoch 

Participating stations 

67 

04:00 

v255q 

2012 March 

ATCA, Cecluna, Hobart, Mopra, Parkes 

77 

04:00 

v255r 

2013 March 

ATCA, Ceduna, HartRAO, Hobart, Parkes 

168 

02:30 

v255s 

2013 June 

ATCA, Ceduna, HartRAO, Hobart, Mopra, Parkes 

226 

18:00 

v255t 

2013 August 

ATCA, Ceduna, HartRAO, Hobart, Mopra, Parkes 

323 

12:00 

v255u 

2013 November 

ATCA, Ceduna, Hobart, Mopra, Parkes 

60 

22:00 

v255v 

2014 March 

ATCA, Ceduna, HartRAO, Hobart, Mopra 


Table 3: Sources in each epoch which were used for ICRF mode delay calibration. 


Epoch 

Source 

name 

RA 

(h m s) 

Dec 

( ° , //) 

Scan 

(min) 

2012 March 

1349-439 

13 52 56.54 

-44 12 40.388 

1:00 

2013 March 

0537-441 

05 38 50.36 

-44 05 08.939 

1:00 

2013 June 

1302-102 

13 05 33.02 

-10 33 19.428 

2:00 

2013 August 

0013-005 

00 16 11.09 

-00 15 12.445 

2:00 

2013 November 

1929+226 

19 31 24.92 

22 43 31.259 

1:15 

2014 March 

0537-441 

05 38 50.36 

-44 05 08.939 

4:00 
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4.2. Phase Referenced Data 

We removed modeled residuals attributed to the ionospheric delay, EOPs and 
parallactic angle before applying the troposphere and clock drift rate corrections which were 
determined from the ICRF observations. The AIPS task CVEL was then used to re-position 
the maser spectrum in the bandpass to account for Doppler shifts due to antenna positions 
and Earth’s rotation specific to each epoch. Figure 1 shows that the G 339.884—1.259, 

6.7 GHz maser spectral peak, has a local standard of rest velocity (ui sr ) of —38.8 km s^ 1 as 
reported by Caswell et al. (2011). 

ACCOR was then used to correct the amplitude of the data for imperfect sampler 
statistics at recording, and also for incorrect amplitude scaling at the correlator which 
was a problem in some versions of the DiFX software used for correlation. Following this, 
we extracted a single autocorrelation scan of the maser spectrum from the most sensitive 
single antenna (Parkes, with the exception of the final epoch) and used ACFIT to scale the 
spectra at all observing stations to this template, based on the nominal system equivalent 
flux density (SEFD) of the antennas 1 . The resultant amplitude gains as a function of time 
were then applied to the maser and quasar datasets. The calibrated peak flux intensity 
for the feature at —35.6 km s -1 in G 339.884—1.259 (the maser component used for phase 
referencing) and J 1706—4600 (the phase referenced quasar used for parallax measurement) 
for each epoch is presented in Table 4. The RMS noise in the quasar image was determined 
from an area of size 1.5 x 1.5 arcsec and from a region of the image where there was no 
emission. The RMS noise for the maser image was determined from an area of the same 
dimensions and from a spectral channel where there was no emission. 

We used a single scan (1 to 2 minutes) on J 1706—4600 for delay calibration to correct 

^ www.atnf.csiro.au/vlbi / documentation / vlbi_antennas / index.html 
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the initial residual clock and instrumental error from correlation. However the quasar 
dataset was not correlated with the same spectral resolution as G 339.884—1.259 (see 
Section 3) and there was a need to account for this by modifying the solutions copied 
to the zoom-band maser data. This crucial step was required in order to prevent the 
introduction of spurious R-L polarization phase differences into the dataset, and the efficacy 
of this procedure (see Appendix B) is demonstrated by the similarity in the phase solutions 
between the different polarizations in Figure 2. Similar phase transfer issues between 
datasets with differing frequency properties have previously been resolved using comparable 
methods (Rioja et al. 2008; Dodson et al. 2014). 

5. Astrometry, Parallax and Proper Motion 

Figure 1 shows that there are multiple strong (>100 Jy) 6.7 GHz methanol maser 
components in G 339.884—1.259 which offer potential spectral channels for astrometry. 

We examined the cross correlation spectra to find the spectral features which showed 
the smallest relative flux density variations across all baselines for the duration of the 
observations, taking this to be an indication that it was an unresolved point source, 
which would enable an accurate position determination. We found the spectral feature at 
—35.6 km s -1 (Figure 3) to be clearly the best choice for G 339.884—1.259. All the other 
peaks exhibited substantial variability across individual baselines suggesting that emission 
is not point-like on milliarcsecond scales, or that there is a blend of emission from different 
locations. We then fringe fit on the maser spectral channel associated with this feature 
before transferring the phase solutions to J 1706—4600. Figure 2 shows a typical plot of 
phase versus time from the 2013 August observations. 

After transferring the phase corrections, we averaged all channels in the quasar dataset 
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Table 4: Measured fluxes and differential fitted positions between the —35.6 km s _1 feature 
in G 339.884—1.259 and J 1706—4600 across all epochs after phase referencing. The formal 
errors in each coordinate are determined from the output of IMFIT. 


J 1706-4600 G 339.884-1.259 

Epoch x Offset Error y Offset Error Flux RMS Flux RMS 

(mas) (mas) (mJy) (Jy) 


2012 March 

4.177 

0.056 

2.998 

0.098 

30.29 

0.55 

413.07 

0.03 

2013 March 

2.667 

0.026 

1.085 

0.023 

28.98 

0.36 

418.17 

0.02 

2013 June 

1.654 

0.003 

0.464 

0.002 

114.84 

0.15 

427.49 

0.02 

2013 August 

1.031 

0.021 

0.208 

0.017 

61.06 

0.65 

462.02 

0.02 

2013 November 

0.727 

0.016 

0.006 

0.012 

62.25 

0.38 

268.95 

0.02 


G339.884-1.259 



20 15 10 5 0 -5 -10 -15 -20 


EW offset (mas) 

Peak flux = 2.6895E+02 JY/BEAM 

Percentage levels = 2.689E-00 * (-90, -64, -32, -8, -4, -2, -1,1,2, 4, 
8,16, 32, 64, 90) 


J1706-4600 



20 15 10 5 0 -5 -10 -15 -20 


EW offset (mas) 

Peak flux = 6.4301 E-02 JY/BEAM 

Percentage levels = 6.430E-04 * (-90, -64, -32, -8, 8,16, 32, 64, 90) 


Fig. 3.— Emission from the phase reference channel corresponding to Vi sr = —35.6 km s _1 in 
G 339.884—1.259 (left) that was strong and showed compact structure. J 1706—4600 (right) 
showed consistent centroid structure dominated by a single peak throughout all epochs. 
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and then imaged the emission using a Gaussian beam of 5.9 x 4.2 mas (average from 
all epochs) (Figure 3). We report detections for J 1706—4600 and J 1654—4812 on VLB1 
baselines, and only the first of these was suitable for astrometry. J 1706—4600 was observed 
in all epochs and appears to be dominated by a single component with no jets. J 1654—4812 
was observed in all epochs except 2012 March and showed variable source structure which 
made it unsuitable for parallax determination. J 1706—4600 has a positional accuracy of 
2.10 mas (Petrov et al. 2011) and shows deviation from point-like structure at levels <10% 
of the peak flux density (Figure 3). J 1654—4812 had an estimated positional accuracy of 
0.4" (McConnell et al. 2012 ) and from our phase referenced images, we are able to present 
updated coordinates for J 1654—4812 to an uncertainty of 2.1 mas in Table 1. We located 
the centroid position of the quasar by fitting a 2D Gaussian to the deconvolved J 1706—4600 
emission. The offset of the emission peak from the centre of the image held was recorded 
for all epochs and we present these data in Table 4. 

The change in the position of the —35.6 km s _1 feature in G 339.884—1.259 with 
respect to J 1706—4600 was modelled independently in right ascension and declination, and 
included the cllipticity of Earth’s orbit. We allowed for systematic sources of uncertainty in 
right ascension and declination in the parallax model and added these in quadrature to the 
formal errors in Table 4. A xt (P er degree of freedom) for the East-West and North-South 
residuals was determined, and we iteratively adjusted the estimated error floors until xl ~ 1 . 
The parallax was measured to be 0.48±0.08 mas corresponding to a distance of 2.1 ^q 3 kpc 
to G 339.884—1.259. The proper motion was found to be fi x = —1.6T0.1 mas y _1 and 
H y = —1.9±0.1 mas y -1 (Figure 4 and Table 5). In order to constrain errors in the measured 
proper motion, we made image cubes of the maser emission for all epochs, and analyzed the 
changes in the distribution from 2012 March to 2014 November. We found that internal 
motions in G 339.884—1.259 can be up to ~0.4 mas y -1 in right ascension and declination. 
These motions dominate over the formal errors in when added in quadrature, 
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and we report the measured proper motion with errors as fi x = — 1.6±0.4 mas y' 1 and 
fiy = —1.9±0.4 mas y _1 . This measured uncertainty corresponds to internal motions of 
~2 km s _1 in the maser emission and is consistent with proper motion estimates of 6.7 GHz 
methanol masers in HMSFRs (e.g. Goddi et al. 2011; Moscadelli & Godcli 2014; Sugiyama 
et ah 2014). A more detailed analysis of the internal motions of the 6.7 GHz emission in 
G 339.884—1.259 is beyond the scope of the current text and will be presented in future 
publications. 

We excluded data from the 2014 March observations in our parallax and proper motion 
analysis as there were significant technical difficulties during that session which prevented 
us from obtaining accurate measurements of the maser-quasar separation for that epoch. 

The LBA has previously been used to measure distances to a number of pulsars at 
1.6 GHz using trigonometric parallax (Dodson et al. 2003; Dellcr et al. 2009b,a). The 
uncertainty in our LBA parallax measurement to G 339.884—1.259 is estimated to be 
80 /iarcsec and is equivalent to the errors of the best LBA southern pulsar parallax 
measurement (Dcller et al. 2009a). The uncertainty from our observations is a factor of 
2 poorer than the parallax of 6.7 GHz methanol masers in ON 1 (0.389 ± 0.045 mas) 
measured by Rygl et al. (2010) using the EVN. The LBA (when operating without 
HartRAO; see Appendix A) and the EVN configuration in Rygl et al. (2010) both have 
a maximum East-West baseline separation of ~2000 km, giving them similar resolution 

Table 5: Parallax distance and proper motion of G 339.884-1.259. 


Maser feature 

(km s' 1 ) 

Distance 

(kpc) 

(mas y' 1 ) 

ID 

(mas y' 1 ) 

-35.6 

o i +0.4 
z ' 1 -0.3 

— 1.6T0.1 

— 1.9T0.1 
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Fig. 4.— Parallax and proper motion of the —35.6 km s _1 reference feature of 
G 339.884—1.259. The expected positions from the fits are indicated with triangular and 
circular markers. Left panel: the sky positions with the first and last epochs labeled. Middle 
panel: East-West (triangles) and North-South (circles) motion of the position offsets and 
best combined parallax and proper motions fits versus time. The models are offset along 
the y-axis for clarity. Right panel: the parallax signature with the best fit proper motions 
removed. 
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capabilities for parallax determination (Reid et al. 2009a). Rygl et al. (2010) reduced 
random errors in their measurements by modelling parallaxes determined from the averaged 
positions of several maser spots. However, given the strong and compact nature of the maser 
spot at —35.6 km s -1 , we are not able to reduce our astrometric errors using this method as 
these will be dominated by systemics from the uncompensated atmosphere. Additionally, 
the the EVN observations of ONI utilise background quasars with separation angles of 1.71° 
and 0.73°, and in comparison J 1706—4600 is separated by 2.48° from G 339.884—1.259 
(Table 1). This would have adversely affected the interpolated phase transfer solutions 
(described in Section 3) and contributed to the larger uncertainty estimate presented here 
(see Appendix A). In considering these factors, we assess that with better atmospheric 
calibration and smaller separation between the maser and quasar, it will be possible to 
attain parallax accuracies of ~20 //areset: using the LBA. 


5.1. Kinematic distance to G 339.884—1.259 

The kinematic distance to sources in the Galactic disk can be determined from a 
model which describes the rotational speed of the disk at the Sun (@ 0 ), distance of the Sun 
from the Galactic centre (Ro), and the measured ui sr . Green & McClure-Grifkths (2011) 
report a kinematic distance to G 339.884—1.259 of 2.6 ± 0.4 kpc using @ 0 = 246 km s -1 , 
R 0 = 8.4 kpc (Reid et al. 2009b) and ui sr = —34.3 km s^ 1 (the mid-point of the 6.7 GHz 
methanol maser emission range Figure 1). Previous studies of the high-mass star formation 
region associated with G 339.884—1.259 (e.g. Ellingsen et al. 1996; De Buizer et al. 2002; 
Dodson 2008) have used a kinematic distance of ~3 kpc to this source. This value was 
determined using earlier models with Galactic rotation speeds of @ 0 — 220 km s^ 1 (see 
Section 6.2 for further comments). 

Using updated Galactic parameters of © 0 = 240 km s -1 , R 0 = 8.34 kpc and solar 
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motion parameters of U Q = 10.70 km s _1 (towards the Galactic centre), V & = 15.60 km s^ 1 
(clockwise and in the direction of Galactic rotation as viewed from the North Galactic Pole) 
and Wq = 8.90 km s -1 (towards the North Galactic Pole) (Reid et al. 2014), we report the 
kinematic distance to the associated CS(2-1) cloud with v\ sr = —31.6 km s -1 (Bronfman 
et al. 1996) to be 2.5±0.5 kpc. This result is comparable to the parallax distance in 
Section 5 but with a large estimated error. Given an uncertainty of ~8 km s -1 in @ 0 (Reid 
et al. 2014), the estimated error in the quoted kinematic distance of 2.5 kpc is doubled to 
~1 kpc. 


5.2. Association of G 339.884—1.259 High-Mass Star Formation Region with 

Scutum Arm 

Our ability to precisely determine the structure of the Milky Way is hampered by our 
location in the midst of its spiral arms. Westerhout (1957); Cohen et al. (1980); Dame et al. 
(2001); Jones et al. (2013); and others have used surveys of Hi and CO molecular clouds 
to identify the Galaxy’s spiral arms from the ensuing longitude-velocity (£ — V) diagrams. 
This method of spiral arm modelling has been used by Xu et al. (2013); Zhang et al. (2013); 
Choi et al. (2014); Reid et al. (2014); Sanna et al. (2014); Sato et al. (2014); Wu et al. 
(2014) to assign HMSFRs to spiral arms, by associating them with molecular clouds in our 
Galaxy. 

Based on the ui sr of the CO spectrum in the region (T. Dame 2014, private 
communication) and the parallax distance in Table 5, we suggest that G 339.884—1.259 
is located in the near edge of the Scutum spiral arm as shown in Figure 5. Sato et al. 
(2014) modeled the Scutum arm based on measurements of 16 HMSFRs, and we show the 
position of G 339.884—1.259 with relation to these in Figure 5. It can be seen that at the 
longitude of G 339.884—1.259, the Sagittarius and Scutum arms are in close proximity, and 
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extrapolating the best information currently available suggests that these two arms may 
merge at lower longitudes. 

Using a log-periodic spiral model, Sato et al. (2014) measured a pitch angle of = 
19.8°±3.1° for the Scutum arm. We have now included G 339.884—1.259 into this model, 
using a source which extends the Galactocentric azimuth by about 10°, giving an updated 
value of = 19.2°±4.1°. 


5.3. Updated coordinates for G 339.884—1.259 and J 1654—4812 

During the observations we used a = 16 h 52 m 04.6700 s , 5 = —46°08 , 34.200" as the 
coordinates for G 339.884—1.259 (Caswell et al. 2011). This is the position of the maser 
spectral peak at —38.7 km s -1 (Figure 1). When applying phase corrections from the 
—35.6 km s _1 feature to J 1706—4600, we can assume that any offset of the quasar image 
from the centre of the held is due to the offset of the phase referenced position with respect 
to the true maser position. We iteratively corrected the maser coordinates in the AIPS 
source table using CLCOR until there was no further improvement in the quasar position 
from the image centre in the 2013 March epoch of observations. It is important to use 
accurate maser coordinates in the data reduction process, to minimize errors in determining 
the quasar position during phase referencing (Reid et al. 2009a). The position corrections 
from the March 2013 epoch of observations were then applied to G 339.884—1.259 for all 
epochs, and we present the updated coordinates corresponding to the —38.7 km s _1 feature 
in Table 1. We find that there is a difference of 0.219" between our measured position 
and that reported by Caswell et al. (2011). The formal errors in the fitted position for 
J 1706—4600 were found to be <0.1 mas (Table 4). As this is an order of magnitude smaller 
than the known positional error of this source (see Section 3), the resultant uncertainty 
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X (kpc) 


Fig. 5.— A face-on view of the Milky Way galaxy showing G 339.884—1.259 in relation to 
other HMSFRs in the Scutum arm (Sato et al. 2014). Sections of arcs of the Perseus (top 
right), Local (top centre with Solar position) and Sagittarius arms are also shown. The 
background circular disks are scaled to approximate the Galactic bar region (~4 kpc) and 
the solar circle (~8 kpc) (see Reid et al. (2014)). 
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in the updated maser position remains as 2.1 mas when error contributions are added in 
quadrature. 

Observations of J 1654—4812 were made using the coordinates a = 16 h 54: m 18.24: s , 

6 = —48 o 13 , 03.7" from McConnell et al. (2012). The updated position of this source when 
measured relative to the corrected G 339.884—1.259 position is also presented in Table 1 to 
an accuracy of 2.1 mas. There is a separation difference of 0.074" between the original and 
updated positions. 


6. Properties of associated high-mass star formation region 

6.1. Peculiar motion 

The proper motions /j x = —1.6T0.4 mas y _1 , /i y = —1.9T0.4 mas y^ 1 and of i>i sr = 
—31.6 km s^ 1 of the CS(2-1) cloud associated with this source (Bronfman et al. 1996), make 
it possible to determine the full 3-dimensional motion of G 339.884—1.259 with respect to 
the Galactic centre. The dynamical model of the Galaxy we use assumes a flat rotation 
curve of the disk with a speed of ©o = 240 km s -1 . This is a reasonable assumption based 
on the recent analysis in Reid et al. (2014). The distance of the Sun to the Galactic centre 
is taken to be R 0 = 8.34 kpc, and assumed to have peculiar motion components U Q = 
10.70 km s“\ Vq = 15.60 km s -1 and W Q = 8.90 km s -1 (Reid et al. 2014). When using 
this model, we find the peculiar motion for G 339.884—1.259 to be U = —4.0±5.9 km s -1 , 
V = 6.47±4.6 km s -1 and W = 10.Oil.2 km s -1 , in a reference frame that is rotating with 
the Galaxy. Hence all components of peculiar velocity are consistent (within estimated 
uncertainty) with the bulk of HMSFRs in the Milky Way. In Reid et al. (2014) there is a 
good fit to the model of spiral arm motions when a RMS of about 5-7 km s^ 1 is assumed 
for each velocity component of HMSFR. This is reasonable for Virial motions of stars in 



giant molecular clouds and so there may not be much evidence for large (>10 km s x ) 
streaming motions in general. 


6.2. Physical constraints on the ionizing star 

High-mass star formation is still not well understood, and there is a vibrant debate 
regarding the processes which result in their formation (e.g. McKee & Tan 2003; Bonnell 
et al. 1998; Garay & Lizano 1999). Accurate distances to star formation regions help to put 
tight constraints on the physical environments from which young high-mass stars are born. 
This includes fundamental attributes such as size, enclosed mass and luminosity. Therefore, 
any accurate model describing the star formation process must also be governed by the 
limits of these physical constraints. 

Previous groups have determined the properties of G 339.884—1.259 based on kinematic 
distance (see Section 5.1). Using values from Ellingsen et al. (1996) and the equations from 
Panagia & Walmsley (1978), we present updated estimates for the electron density (n e ), 
mass of ionized hydrogen (Mhh), excitation parameter ( U ) and the Lyman continuum 
photon flux (Nl) in Table 6. Based on log Nl = 45.6 s _1 in Table 6, we find that 
continuum emission is consistent for a core object to be a B1 type star (Panagia 1973). 
This classification implies that G 339.884—1.259 is relatively small for a young high-mass 
star. It has been proposed that the luminosity of 6.7 GHz methanol masers increases 
as the associated young stellar object evolves (e.g. Breen et al. 2010). In this scenario 
sources such as G 339.884—1.259, which is amongst the most luminous 6.7 GHz methanol 
masers in the galaxy and is associated with a number of rare class II methanol transitions 
(Krishnan et al. 2013; Ellingsen et al. 2013), are close to the end of the methanol-maser 
phase of high-mass star formation. Recently, Urquhart et al. (2015) have put forward an 
alternative interpretation, that the methanol maser luminosity has a closer dependence on 
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the bolometric luminosity of the associated high-mass star than its evolutionary phase. 

In this case, we would expect G 339.884—1.259 to be associated with a high-mass young 
stellar object towards the upper end of the mass-range. However, this is not the case and 
demonstrates that high luminosity, multiple transition methanol maser emission need not 
be associated with the most massive O-type young stars. A single example is insufficient 
to resolve the issue of whether evolutionary stage or stellar mass plays the primary role in 
determining class II methanol maser luminosity, however, G 339.884—1.259 shows greater 
consistency with the expectations of the Breen et al. (2010) evolutionary hypothesis. 


7. Conclusion 

We are currently conducting a large project using the LBA to measure the positions 
of thirty, 6.7 GHz methanol masers relative to background quasars, in the Southern 
hemisphere to sub-milliarcsecond accuracy. These measurements will be used to determine 
their distances using trigonometric parallax. The source list contains many prominent 
HMSFRs in the third and fourth quadrants of the Milky Way galaxy, where the LBA’s 
performance is unmatched in its astrometric capabilities as a VLBI instrument. In this 
paper, we have shown the potential of this project by successfully making the first parallax 
measurements to a southern 6.7 GHz methanol maser source. The parallax of 0.48±0.08 mas 

Table 6: Physical parameters of G 339.884—1.259 as described in Section 6.2 and adjusted 
to a preferred distance of 2.1 kpc. 


n e 

Mh„ 

U 

log N l 

Spectral 

(cm -3 ) 

(M 0 ) 

(pc cm -2 ) 

(s -1 ) 

type 

3.1 x 10 4 

6 x 10 -4 

5.7 

45.6 

B1 
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to G 339.884—1.259 corresponds to a distance of 2.1^0 3 kpc. In combining this result with 
measurements of other HMSFRs in Sato et al. (2014), we place G 339.884—1.259 at the near 
edge of the Scutum spiral arm of the Milky Way, and determine an updated pitch angle of 
-0 = 19.2°±4.1° for this arm. 

We have used the parallax distance to update the estimated physical parameters for the 
G339.884—1.259 star formation region and now classify it to be of spectral type Bl. The 
young stellar object associated with G 339.884—1.259 is relatively small for a “high-mass” 
object and is a strong example of the uncorrelated nature between stellar mass and intensity 
of 6.7 GHz maser emission. 

Results from VLBI astrometric observations including the BeSSeL survey are continuing 
to shape the way we view our Galaxy, by clearly revealing its spiral arm structure, rotation 
dynamics and mass. These results are based on observations which are concentrated in the 
first and second quadrants of the Milky Way (Reid et al. 2014), and the inclusion of results 
from the LBA is vital to ensure that sources are sampled from across a broad range of 
Galactic longitudes in order for a complete picture of the Galaxy to emerge. 
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A. Tropospheric and Ionospheric Contributions to Astrometric Accuracy 

We found that the parallax measurement with the smallest error floors was obtained 
when troposphere corrections were applied to the 2013 March and 2013 August epochs 
only. When the corrections were incorporated for all epochs, the parallax was measured 
to be 0.58±0.11 mas. When we excluded the tropospheric corrections for all epochs the 
measured trigonometric parallax was 0.50±0.13 mas. We attribute the effectiveness of 
applying this calibration in each epoch to limitations that we encountered in estimating 
both the tropospheric and ionospheric delay. 


A.l. Troposphere Calibration 

The multiband delay solutions obtained from 4 x 16 MHz IFs were found to be 
significantly noisier than those obtained when only 2 x 16 MHz IFs were used. As a 
result, the multiband delays for our observations were determined from a total bandwidth 
separation of only 32 MHz instead of the intended 400 MHz. This led to a loss in the 
accuracy of the delay measurement and affected our ability to model the troposphere path 
length contribution. 

Deller et al. (2009a) demonstrated that the available TEC models, for ionosphere delay 
correction in the Southern hemisphere, are sometimes in error due to the lower density of 
GPS stations in this region. The ionospheric delay is dispersive whereas the tropospheric 
delay is not. Hence uncorrected ionospheric delays which are inadvertently handled as 
non-dispersive tropospheric delays can result in erroneous phase corrections (Rygl et al. 
2010). This could be the reason why the search for multiband delays over the highly 
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sensitive 400 MHz range was noisy, and we detail the ionospheric path length contribution 
in the next subsection. 

In addition to the troposphere path length, the clock drift over the course of the 
observations was also resolved from the multiband delays. The clock model assumes that 
the H-maser at each station had a delay drift which varied linearly with time. Figure Al 
shows the residual differences between the measured clock and modeled offsets for the 
duration of an observation. The squares, circles and crosses represent (respectively) the 
data, the model and residuals between the two. The near zero scatter of the residuals 
indicates that a linear model for the clock drift at each station was successful. The RMS 
noise of the residuals were found to be between 0.03 to 0.3 nanoseconds across all baselines 
and epochs. 

The hetrogeneous nature of the array means that some antennas were able to 
participate in a relatively small number of ICRF scans, Parkes due to a combination of 
slow slew rates, limited elevation coverage and HartRAO due to its distance from the 
rest of the LBA antennas. We were not able to include HartRAO in the majority of the 
ICRF observations as the quasars which had risen at the Australian stations were often set 
here. As a result, we were unable to determine the clock drift for HartRAO and have not 
included observations from this station to determine the results given in Table 5. We are 
currently investigating alternative methods to determine the clock drift rate from HartRAO 
for inclusion in the future. 


A.2. Ionosphere calibration 

Given the observing parameters (6.7 GHz, 4 minute cycle time, 2.48° separation 
between the calibrator and the target, 3 cm residual zenith path length (Reid et al. 1999) 
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Fig. Al.— A sample of multi-band delay solutions from Ceduna (CD), Hobart (HO), Mopra 
(MP) and Parkes (PA) to the ATCA (AT), determined from the ICRF mode observations 
in 2013 June. The squares, circles and crosses represent (respectively) the data, the model 
and residuals between the two. 
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and residual ionospheric content of 6 TECU (Ho et al. 1997)) we can estimate the expected 
error contributions from static and dynamic components of the troposphere and the 
ionosphere using the formulae in Asaki et al. (2007). This predicts a dynamic tropospheric 
phase error of 28°, a static tropospheric phase error of 15°, a dynamic ionospheric phase 
error of 5° and a static ionospheric phase error of 21°. The geodetic blocks typically reduce 
the residual zenith path length to about 1 cm, which would decrease the static tropospheric 
phase error to 5°. Dynamic errors, because of the short timescale on which these operate, 
reduce the measured flux density of the targeted source without having a large effect on 
the positional centroid. This clearly leaves the residual static ionosphere as the largest 
uncorrected effect. Static residuals will introduce a shift in the observed position, as 
measured on any single baseline. For arrays with 10 or so antennas these shifts average out 
somewhat, but the LB A with 5 to 6 antennas is more vulnerable to this effect. 

Tuning the observational parameters, such as the cycle time and the separation between 
sources, is the most effective method to reduce these error contributions. For example, 
reducing the cycle time to 60 seconds would diminish the dynamic phase errors to less than 
10° and decreasing the separation between the calibrator and the target to 1° would cut 
back the static phase errors to less than 10°. However, the array sensitivity places limits the 
minimum useful scan duration and the separation of the closest suitable phase calibrator. 
In this case we were not able to have shorter scans or closer calibrators. 

For observations at higher frequencies, such as 22 GHz, the contribution of the residual 
zenith path length dominates. Therefore the geodetic blocks are absolutely essential for the 
measurements in Reid et al. (2014); Honma et al. (2012) as they reduce the typical errors 
from 3 cm to about 1 cm. However at 6.7 GHz the residual ionospheric contribution of 6 
TECU is equivalent to 5.3 cm, as opposed to 0.5 cm at 22 GHz, and there is currently no 
established strategy for minimizing these contributions. Alternative approaches to lower the 
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residual ionospheric contribution are currently being investigated and will be the subject of 
future publications. 


B. Delay Calibration for Phase Referencing Data 

We used FRING to determine the visibility phases and group delay from an individual 
scan on the delay calibrator J 1706—4600 taken from the continuum mode data. It is 
assumed that the delay r calculated by FRING is a constant for the IF and can be used to 
determine a phase correction A <f> for a frequency channel with width Av and with respect 
to the lower-band edge given by 


A <p = tAv 


The 2 MHz maser zoom-band is a sub-band of the IF used to determine the manual 
delay, but it has a different lower band-edge (offset by R 0 //) to the continuum data and a 
different spectral channel width bv. To apply the appropriate phase correction to the maser 
reference channel it was necessary to edit the FQ table in AIPS so that the spectral channel 
width bv' times the channel number a pr corresponds to the frequency of the maser phase 
reference channel v pr , referenced to the lower band-edge of the continuum mode data given 
by 
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It was also necessary to edit the total bandwidth parameter in the FQ table for the 
maser data so that it matched the bandwidth of the continuum mode data. 

This step can be avoided if the delay calibrator scans are correlated for both continuum 
and zoom mode configurations and it is recommended that this be the general procedure 
in future LBA observations. Equivalent solutions for A0 and r can be obtained for both 
datasets if an identical time range is used in FRING. An alternative approach would be 
to correlate the continuum data into several channels such that one of these corresponds 
exactly to the band coverage of the zoom mode data. The multi-band group delay solutions 
obtained from the continuum mode data could then be used to find the phase solution in the 
channel corresponding to the zoom band. This solution could then be directly transferred 
to the zoom band data. 
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